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ABSTRACT 



Context. We present medium spectral resolution near-infrared (NIR) HK-hand spectra for eight low redshift (z < 0.06) radio galaxies 
to study the NIR stellar properties of their host galaxies. The sample was selected from the radio galaxy sample imaged previously in 
the B- and R-hand. They were found to be bluer than inactive elliptical galaxies, possibly indicating a recent star formation episode. 
As a homogeneous comparison sample, we used nine inactive elliptical galaxies that were observed with the same telescope and 
detector with similar resolution and wavelength range. 

Aims. The aim of the study is, by using the advantage of NIR absorption features, to compare the NIR spectral properties of radio 
galaxies to those of inactive early-type galaxies and, furthermore, produce the first NIR HK-hand spectra for low redshift radio galax- 
ies. 

Methods. For the radio galaxy and inactive elliptical samples, spectral indices of several diagnostic absorption features, namely 
Sil(1.589yum), CO(1.619/im) in the //-band and Nal(2.207yum), Cal(2.263yum), CO(>2.29//m) in the /f-band, were measured. The 
strength of absorption lines depends on the luminosity and/or temperature of stars and, therefore, spectral indices can be used to trace 
the stellar population of galaxies. To characterize the age of the populations, the measured EWs of the absorption features were fitted 
with the corresponding theoretical evolutionary curves of the EWs calculated by the stellar synthesis model. 

Results. On average, EW(CO 2.29) of radio galaxies is somewhat greater than that of inactive ellipticals. Most likely, EW(CO 2.29) 
is not significantly affected by dilution, and thus indicating that elliptical galaxies containing AGN are in a different stage in their 
evolution than inactive ellipticals. This is also supported by comparing other NIR absorption line features, such as Cal and Nal, with 
each other. Based on the theoretical evolutionary curves of EWs absorption features are consistent with the intermediate age stellar 
population, suggesting that host galaxies contain both an old and intermediate age components. Intermediate age population is also 
consistent with previous optical spectroscopy studies, which have shown evidence on the intermediate age (~2 Gyr) stellar population 
of radio galaxies, and also in some of the early-type galaxies. 

Conclusions. Intermediate stellar population component indicates that radio galaxies have experienced a star formation epoch rel- 
atively recently. The existence of intermediate stellar population is a link between the star formation episode, possibly induced by 
interaction or merging event, and the triggering of the nuclear activity. 
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1. Introduction 
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Recent observations have shown that most bulge-dominated 
galaxies con tain a supermassive black hole in their center (e.g., 
ICombesll2005 ). This indicates that nuclear activity phenomenon 
can be a transient phase in evolution history of a galaxy that 
can be understood in terms of the central black hole fueled 
by in-falling gas toward the central region. The in-falling gas 
might be triggered either by merging processes of the spheroid 
formation at high redshi ft or during milder interac tion of the 
galaxies at lower redshift (I Cavaliere & Vittorinil2000l) . Dynamic 
interaction of an early-type galaxy with a gas rich late-type 
galaxy may often induce strong star formation as indicated by 
spectroscopic studies of local field early -type galaxies belong- 
ing to pairs or exhibi ting shells (e.g., iLonghetti et al.l l2000t 
iTantalo & Chi osi"2004') and theoretical studies of the merging 
process (e.g., Koiima & Noguchi 1997). Thus, the same mech- 



anism may be responsible for both the nuclear activity and star 
formation. Characterizing the stellar populations in the central 
regions of galaxies, where supermassive black holes and recent 
star formation may co-exist, provides essential information for 
probing their origin and evolution. AGN feedback provides a 
mechanism to quench star formation on a short timescale (e.g., 
Springel et al. 2005) and to re plicate the observed b imodal color 
distribution of galaxies (e.g.. IStrateva et al]|2001h . The energy 
released by the AGN expels enough gas to quench both star for- 
mation and further black hole growth and, thus, determines the 
duration of starburst and AGN phase. 

Most luminous AGN are hos ted by large and luminous 
(massive ) ell i ptical g alaxies (e.g . , iFalomp & Kotilaineri '199Jt 



Govoni et~an 120001; iHeidt et ali 120041; 



_ - - - . : ■ iK otilainen & Falomo| 

l2004c iKotilainen et al.1 120051; iHvvonen et alj |2007b). There is 
increasing evidence that the host galaxies of vario us types of 
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AGN are pi t ier tn a n inactive ellipticals ( uovoni et ai 
ISchade et all 120001; lOrndahl et al] 120031; iSanchez et al. 



bluer than inactive ellipticals ( Govoni et al 



2OO0; 



2004 
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JahnkeetalJ 12004 iKotilainen & Falomol I200I iHvvonen etal] 



2007al) and they do not follow the color-magnitude relation of in- 



active ellipticals (Peletier et al. 1990) indicating a recent star for- 
mation episode. This scenario is also supported by optical spec- 
troscopic studies of nea rby radio galaxies (RG) dRaimann et al.l 
120051: IHoU et al.l l2007h . It is worth to note that the observed 
blue colors are not affected by the dust content of the galaxy, 
because in the presence of dust reddening the intrinsic color 
of the galaxy would be even bluer Inactive elliptical galaxies 
have experienced a relatively short sta r formation episode a t 
z ~2 and evolved passively ever since jCarretero et al.l l2"007h . 
However, the distribution of the ages of stellar population of in- 
active ellipticals is inhomogeneous, often showing an interme- 
diate age population (~3 Gyr), in addition to the energetically 
dominant, very o ld (~15 Gyr) populations dYamada et al.ll2006t 
ISilva et aril2008h . It is important to note that there is a down- 
sizing effect, in the sense that the intermediate age populations 
are mostly found in early-type galaxies with low velocity dis- 
persion, i.e. low mass, while the most massive galaxies are all 
very old. However, there are also indications that galaxies in 
den ser enviroments (clusters) are slightly older than field galax- 
ies jTrager et alj|2000l: iKuntschner et al.ll2002l) . It is likely that 
the old stellar population is also dominating the stellar content 
of AGN hosts which has undergone a single recent star forma- 
tion burst with only a small fraction of mass involved in the 
episode. The broad distribution of the colors most probably indi- 
cates an object-to-object difference in the age since the last star 
formation episode. The existence and the age of the young stel- 
lar population is a link between the interaction and/or merging 
event and the star formation episode and the triggering of the 
nuclear activity. Unfortunately, age estimates are not straightfor- 
ward because of the metallicity effect. Strength of CO absorption 
lines of old and metal-rich stellar populations is similar to young 
starbursts with lower metallicity. Also, the cooler the younger - 
relation does not always hold because temperature depends on 
metallicity. 

Luminosity of the majority of AGN is dominated by bright 
nuclear source contaminating the luminosity of the underlying 
host galaxy. Thus, spectroscopic observations of the stellar con- 
tent of the host is diffi c ult. Such studies ( e .g., iBoroson et alj 
1985"; 'Oliva et al.l [l99l [Origlia et al.1 fT997t lOliva et all [1991 
me^ 



Mobasher & James] 120001) have usually been limited by small 
heterogeneous samples, arrays and wavelength range. These 
studies hav e only recen t ly been follow e d by s ystematic stud- 
ies by e.g. [Nolan et alj (|2001[) . iJahnkd ( l2002h . iRaimann et all 
(l2005 h and ^ Silva et al.l (l2008h . According to the unified model 
of radio-loud AGN, RGs are seen edge-on so that the bright 
central region is obscured making spectroscopy of the inner 
part of the host galaxy much easier UV/optical spectroscopy 
of the inner region of RGs suggests that a large number of 
RGs contain a significant contribution from a young stellar 
population (e.g., lAretxaga et al. 2001; Tadhunter et al. 2002; 
IWills et al.1 12001 Ikaimann et al J 120051: IHoltet al.ll2007i) . prob- 
ably related to star formation episode in the near past. Optical 
Lick spectral indices have usually been used in stellar popula- 
tion studies dWorthev & Ottavianilll997l) . but they have prob- 
lems due to blended features and various population of galax- 
ies. Furthermore, in the optical wavelengths, spectral features 
from the young stellar population are difficult to observe due to 
possible presence of other emission components, such as neb- 
ular continuum, and direct and scattered AGN emission. NIR 
spectroscopy has advantages compared to optical, because RGB 
stars dominate at 2jum and the wavelength region contains many 
diagnostic stellar absorption lines (e.g., Sil 1.589//m, CO(6-3) 



L619//m, Nal 2.207jum, Cal 2.263//m and CO(2-0) bandhead 
>2.29jum) which can be used as indicators of stellar population 
in terms of their temperature and/or luminosity. 

NIR spectroscopy has mostly been used to study high red- 
shift RGs and t here are only few stu dies for individual so urces 
at low redshift dBellamv et a l. 2003; Tadhunter et"ani2003l) . We 
present NIR spectroscopy for eight nearby (z < 0.06) RGs (one 
of them, PKS 0521-36 is also classified as BL Lac object) with 
reasonable spatial resolution (~ Ikpc) and medium spectral res- 
olution {R ~ 1000) to study the NIR properties and ages of their 
stellar populations compared with the sample of nine inactive 
ellipticals. Except for individual sources (e.g., iTadhunter et al.l 
,2002) these are the first NIR spectra of a sizable sample of low 
redshift R Gs. The RG sa n iple w as previously observed in R- and 
B-band by 'Fasan o et alj (Il996h . and detailed B-R color infor- 
mation are availa ble for them. Th e sample was selected among 
the objects of iGovoni et al.l (1200 0) that are bluer than inactive 
elliptical galaxies and do not follow the B-R color-magnitude 
relation of inactive ellipticals. Long-slit spectra of RGs are ex- 
tracted across the nucleus to detect all the important diagnostic 
absorption lines, and extended emission of the host galaxy can 
be studied to large distance from the nucleus. Our aim is to study 
the stellar population of RGs based on the NIR stellar absorption 
indices, and to compare the indices of RGs with inactive coun- 
terparts of similar morphology types. This will have important 
implications for understanding nuclear star formation histories 
and the evolution of RGs. 

In Section 2, we describe the sample, observations, data re- 
duction, and methods of analysis. In Section 3, we present the 
results and discussion concerning the properties of the galaxies. 
Summary and conclusions are given in Section 4. 



2. Observations, data reduction, and analysis 

The observations of the RGs were carried out in September 2005 
at the ESO New Technology Telescope (NTT) using SOFI cam- 
era with pixel scale 0.288 arcsec pix"'. The medium resolution 
{R ~ 1000) long-slit ///iT-band spectra were taken across the nu- 
cleus with red grism and slit width 1 .0 arcsec. The wavelength 
range of the spectra is from 1.5 to 2.5jum with typical integration 
time of 1440 s. The spectra of each target were taken with pairs 
using ABBA observing cycle and moving the target in the slit. 

Data reduction was performed using IRAlfl Pairs were sub- 
tracted from each other to eliminate background sky emission 
and then divided by the flat-field image. Bad pixels and cosmic 
rays were masked out. Wavelength calibration was done using 
Xe-arc lamp calibration frames. Galaxy spectra were flux cali- 
brated and divided by spectrophotometric standard star and then 
averaged. The nuclear (1.5 arcsec corresponding to ~1.2 kpc in 
physical scale) and off-nuclear extended emission spectra were 
extracted for each target. In order to achieve sufficient S/N-ratio, 
the off'-nuclear spectra were added together from both sides of 
the nucleus (1.5 arcsec/~1.2kpc at each side). The //-band spec- 
tra were normalized by fitting the continuum level to both sides 
of the absorption features. In the /T-band the continuum was fit- 
ted only to the shortward wavelengths from the CO bandhead. 
The journal of the RG sample is presented in Table[T] 

To have a reliable comparison between RGs and inactive 
ellipticals it is important that both samples have similar spec- 



' IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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tral resolution and aperture. As an inactive comparison sample 
we used the sample obtained and kindly provided for us by V. 
Ivanov (see also Cesetti et al. (2008)). The comparison sample 
consists of 7 targets that were also observed by NTT/SOFI in 
March/April 2006 with an identical instrument configurations. 
These galaxies are classified as giant ellipticals with morphol- 
ogy type r = -5 but, in order to have morphology type T 3 

galaxies also in the comparison sample, we included two galax- 
ies from the inactive early-type galaxy survey observed previ- 
ously by us with NTT/SOFI using the same set of parameters 
and the reduction process as for the RG sample (Reunanen et 
al., in prep.). This data altogether gives us optimal spectral res- 
olution matched sample for stellar population comparison. The 
journal of the inactive elliptical galaxy sample is presented in 
Table 121 

The equivalent widths (EW) of the spectral lines were de- 
rived by integrating over the lines of the normalized spectra. 
The integration limits for the particular line were the same for 
all targets. The accuracy of the EWs is limited by the difficulty 
of determination of the continuum level in the vicinity of the par- 
ticular line. Especially this holds for the //-band which is rich in 
spectral features, and for the CO(2-0) bandhead in the /T-band 
because the continuum level has to be extrapolated longward of 
2.295yum. Furthermore, it is worth to keep in mind, that the stel- 
lar features are blended due to atomic and molecular features 
ace & Hinkle 1996). There is no velocity dispersion data 
available for the RG sample, and we have assumed that they are 
roughly in agreement with those of inactive ellipticals used in 
this work. However, it is worth to note that the measured EWs 
of the indices are evidently upper limits due to spreading caused 
by the velocity dispersion. The possible dependence of the EWs 
on velocity dispersion and resolution was c hecked by degrad- 
ing the medium resolution stellar spectra from lWallace & Hinkld 
(Il997h to the lower resolution applicable for the RG sample, 
and measuring the indices of all the five diagnostic features for 
both resolutions. This comparison showed that the indices do 
not significantly depend on the degradation, in the sense that 
all the measured indices agree within the errors. In the normal- 
ized spectrum, a straight line was fit to the continuum points 
around the spectral feature. Central wavelength ranges for each 
feature were 1.585-1.593//mfor Sil, 1.616-1.627//mforCO(6-3) 
and 2.2924-2. 2977yum for CO(2-0). For a direct co mparison, the 
same w ave length ranges wer e adopted as those in ' Origlia et al.l 
( Il993h and lSilva et al] (|2008|) . The errors of the EWs are domi- 
nated by the uncertainties in fitting the continuum level and are 
typically ~0.5 A in //-band and /T-band metallic lines, and ~ 1 .OA 
for the /T-band CO(2-0) feature. 

3. Results 

The ///T-band off-nuclear and nuclear spectra for RG sample are 
presented in Fig.[T] In all off-nuclear and nuclear spectra of RGs, 
//-band spectral absorption features from stellar population (Sil, 
CO(6-3)) are clearly visible as are in the most cases the fea- 
tures Nal, Cal, CO(2-0) in the /T-band. The EWs of all diagnos- 
tic spectral features can be measured for all inactive elliptical 
galaxies, except Sil in two cases (NGC 4472, NGC 4564). The 
EWs of the absorption lines for RGs and inactive ellipticals are 
presented in Table [3] 

By combining each individual redshift corrected spectra to- 
gether, we also created a representative average spectrum for 
RGs and inactive ellipticals (PKS 0521-36 was not included in 
the composite spectrum of RGs) which are shown in Fig. |2] 
These spectra can find applications as template spectra for the 



Table 1. Journal of the radio galaxy sample." 



Target 




T 


B-R 


Mk 




Date 










mag 


min 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


PKS 0023-33 


0.0498 


-5 


1.47 


11.3 


32 


27/09/2005 


3C 29 


0.0450 




1.57 


11.0 


24 


26/09/2005 


NGC 612 


0.0298 


-1.2 


1.59 


10.1 


24 


25/09/2005 


ESQ 552-G14 


0.0317 


-2.7 


1.56 


10.6 


24 


27/09/2005 


PKS 0521-36 


0.0553 




1.48 


11.0 


32 


25/09/2005 


ESQ 528-G36 


0.0408 


-5 


1.48 


10.8 


24 


27/09/2005 


NGC 6998 


0.0397 


-5 


1.84 


10.8 


32 


26/09/2005 


PKS 2158-380 


0.0333 




1.13 


114 


32 


27/09/2005 



" Column (1) gives the name of the target; (2) redshift; (3) morphol- 
ogy type; (4) the S - colour l lFasano et al]|l996h ; (5) /iT-band 
magnitude from 2MASS in 5 arcsec radius aperture; (6) integration 
time, and (7) the date of the observation. 



Table 2. Journal of the elliptical galaxy sample.'' 



Target 


z 


T 


T 

exp 


Date 


(1) 


(2) 


(3) (4) 


(5) 


NGC 2640 


0.00350 


-3 


'il 


09/03/2006 


NGC 4478 


0.00469 


-5 


40 


17/04/2006 


NGC 4472 


0.00333 


-5 


33 


17/04/2006 


NGC 4546 


0.00350 


-3 


32 


11/03/2006 


NGC 4564 


0.00377 


-5 


30 


17/04/2006 


NGC 4621 


0.00143 


-5 


20 


18/04/2006 


NGC 4649 


0.00365 


-5 


33 


17/04/2006 


NGC 4697 


0.00415 


-5 


23 


17/04/2006 


NGC 5576 


0.00505 


-5 


25 


17/04/2006 



Column (1) gives the name of the target; (2) redshift; (3) morphol- 
ogy type; (4) integration time, and (5) the date of the observation. 



underlying galaxy which can be subtracted from the spectra of 
composite stellar systems (i.e. galaxies hosting an AGN). 

3. 7 . Near-infrared line indices as a stellar population 
indicators 

NIR spectra can be used as a tool for studying stellar popula- 
tion in galaxies and star clusters, because the spectral domain 
contains many strong stellar absorpti on features . Especially, the 
/T-band has been widely s tudied (e.g..iMerril & Ridgwav 19791: 
iKleinmann & Hailll986t [Wallace & Hinkld 1 1997h . and it is a 
natural choic e for studying redd ened (du sty) galaxies at low 
redshift (e.g., iMobasher & Jame s 2000; M annucci etani200lh . 
However, //-band is even more suitable for population stud- 
ies because in this band non-stellar emission, mostly caused by 
dust, has smaller effect than in the /T-band but //-band domain 
has been studied systematically only relatively recently (e.g., 
Origlia et al. 1993; Oliva et al. 1995; Fremaux et a n i2007h . 

The most evident absorption features in the /T-band are 
the molecular CO(2-0) bandhead (>2.295/im), atomic Br-y 
(2.166yum), Nal (2.207/zm), and Cal (2.263yum) lines. Br-y is 
the strongest absorption feature within stars earlier than K5, but 
becomes undetectable in stars later than K5 ( Kleinmann & Hall 
1986), while Nal is the stronges t atomic feature when effective 
temperature is Te/f < 3400/: (lAli et al.lll99l) . //-band has a 
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Fig. 1. The HK-hand ofF-nuclear (upper) and nuclear (lower) spectra of radio galaxies. The most prominent diagnostic spectral 
absorption features together with supernova activity indicator emission lines (H2 and [Fell]) are marked by dotted vertical lines. 
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Fig. 2. The HK-hand composite spectra for radio galaxy (upper panels) and inactive elliptical galaxy (lower panels) samples. The 
most prominent diagnostic spectral features are marked by dotted vertical hnes. 



Hyvonen et al.: The stellar population of low redshift radio galaxies 



5 



Table 3. EWs of the absorption lines for off-nuclear and nuclear 
spectra of radio galaxies and inactive elliptical galaxies in units 
of A. 



Target 
A [jum] 



Sil CO(6-3) Nal Cal CO(2-0) 
1.589 1.619 2.207 2.263 > 2.295 



RGs: off-nuclear 
PKS 0023-33 
3C 29 
NGC 612 
ESQ 552-G14 
PKS 0521-36 
ESQ 528-G36 
NGC 6998 
PKS 2158-380 

RGs: nuclear 
PKS 0023-33 
3C 29 
NGC 612 
ESQ 552-G14 
PKS 0521-36 
ESQ 528-G36 
NGC 6998 
PKS 2158-380 

Ellipticals: 
NGC 2640 
NGC 4478 
NGC 4472 
NGC 4546 
NGC 4564 
NGC 4621 
NGC 4649 
NGC 4697 
NGC 5576 



-1.2 
-2.7 

-5 
-5 



-5 

-1.2 
-2.7 



5.3 
4.5 
2.6 
3.4 
1.2 
2.3 
2.8 
4.8 



3.1 
2.5 
2.6 
2.3 
0.8 
1.8 
2.3 
3.8 



5.1 
1.7 

4.8 

2.5 
2.1 
2.1 
2.2 



5.5 
4.6 
4.5 
4.5 
2.5 
3.3 
4.2 
3.5 



6.0 
3.8 
3.8 
5.1 
3.1 
5.8 
3.3 
5.1 



3.1 
5.0 
4.5 
4.5 
4.0 
6.5 
5.1 
5.5 
5.9 



4.2 
3.4 
4.4 
1.1 
4.4 
3.9 
4.9 



4.5 

3.9 
2.5 

3.7 
5.3 
3.9 



5.1 
3.5 
5.5 
4.5 
5.5 
5.2 
4.9 
4.1 
4.0 



3.7 
3.5 
3.0 
2.9 
2.5 
3.1 
3.8 
3.3 



2.9 

3.4 
2.3 

7.7 
5.8 
3.4 



3.1 
3.7 
2.7 
2.5 
2.8 
5.6 
2.7 
3.5 
2.4 



17.1 

20.5 

18.8 

15.1 

5.4 

20.7 

29.0 

20.6 



17.5 

20.1 
14.9 

13.3 
16.5 
14.6 



20.7 
12.5 
17.4 
17.3 
17.4 
16.6 
17.2 
16.3 
14.2 



very complex absorption line pattern due to a number of metal- 
lic and molecular lines, especially in cool stars. The most studied 
//-band features ar e atomic Sil (1.589 jum) and molecular CO(6- 
3) (1 .619yum) (e.g.. lOriglia et al.lll993l) . It is well known, that the 
Nal, Cal and CO feature s in K and M giant sta r s become stronger 
with red der J-K color dRamirez et al.l 119971: iForster Schreibeii 
[2000 : Frogel et al.ll2001h. At a given J-K, giants i n more metal 
rich c lusters have stronger Nal and Cal features (iFrogel et al.l 
1200 lb . In the NIR, the integrated light is dominated by stars that 
are very different from those that contribute significantly at visi- 
ble wavelengths. The first overtone '^CO bandheads are the most 
prominent features in the NIR spectra of old and intermediate 
age populations, and so are potentially prime diagnostics of stel- 
lar content. The depth of these features are predicted to be sen- 
sitive to population parameters (e.g., Vazquez et al. 2003). Nal 
and Cal indices sample transitions from more than one element, 
complicating the sensitivity of these lines to chemical abun- 
dance. Such contamination is also common among optical line 
indices fe.g.. lWorthev et al]|1994l) . Whereas the integrated light 
at visible wavelengths comes from a mix of stellar types, the K- 
band light from all but the youngest stellar systems is dominated 
by evolved giants and AGB (Asymptotic Giant Branch) stars, 
whose photometric and spectroscopic properties are sensitive to 
age and metalUcity (iFrogel et al.lll978l:lMarastonl200 5). The rel- 
ative strengths of the NIR features can be used to compare the 
stellar contents of the central regions of galaxies, and to provide 
insights into their history. 



Currently there are no self-consistent theoretical spectral 
synthesis models for the interpretation of NIR spectra of galax- 
ies. Thus, one has to rely on e.g. their com parison with high- 
resolution NIR stellar spectral atlases (e.g.. IWallace & Hinkld 
1996) . The usefulness of stellar absorption lines as tools for trac- 
ing the stellar population of galaxies is based on the fact that 
particular lines are sensitive to the spectral type (temperature) 
and/or lu minosity ( gravity) of the st ar (Kleinmann & Hall 1983 
lAli et al.l ll995: Ra mirez et al.lll997h . CO hues decrease with in- 
creasing temperature due to variation of the CO/C abundance 
ratio, being detectable at temperatures lower than T^/f ~ 3500K 
and increasing with luminosity. Thus, /T-band CO bandhead is 
very strong in young giant and supergiant stars (10 Myr-100 
Myr) and strong in cool asymptot ic giant branch (AGB) stars 



iv mptot ic 

(100 Myr-1 Gvr:' 01iva et all (1 19951) ). while it is weaker in older 
population. This make s CO bandhea d suitable to trace recent star 
formation in galaxies (Mavval ll997l) . However, some line ratios, 
such as EW(CO 1.62)/EW(SiI 1.59) and EW(CO 1.62)/EW(CO 
2.29), are even bet ter temperature indicators than single lines 
jOriglia et alJll993l) . Especially, the former ratio is useful be- 
cause the two lines have very similar wavelengths and therefore 
the ratio is not sensitive to dilution of the emission of the hot gas. 



3.2. Near-infrared line indices of radio galaxies 

Non-stellar thermal, such as dust surrounding young star form- 
ing regions, and/or non-thermal dilution reduce the intrinsic val- 
ues of EWs originating from the stellar population. However, 
the dust observed in NIR is close to its sublimation temperature, 
and in such cases HII emission line series are expected to be ob- 
served. The comparison of absorption indices in different filters 
depends on the level of dilution by the non-stellar component 
that can be different in H- and /iT-band. The dilution fraction of 
the continuum emission at 1 .6/^m can be estimated from the plot 
of EW(CO 1.62) vs. EW(CO 1.62)/EW(Si 1.59) which is shown 
in Fig. [3] Objects with diluted stellar features lie below the lo- 
cus occupied by stars, because they have a shallower CO(1.62) 
index, while EW(CO 1.62)/EW(Si 1.59) is not significantly af- 
fected by dilution since two the features lie close in wavelength. 
The fraction of non-stellar continuum is simply given by vertical 
displacement of the point in the diagram. The dilution at 2.3yL(m 
can be determined in a similar way from the plot of EW(CO 
1.62) vs. EW(CO 1.62)/EW(CO 2.29) shown in Fig.|4] In that 
case, objects with significant non-stellar continua lie on the right 
hand side from the enclosed region. Note, that one first needs to 
correct EW(CO 1.62) for dilution before determining the non- 
stellar fraction from the horizontal displacement in the diagram. 
In the H- and /T-band, all off-nuclear RGs have either no or only 
slight dilution, while the BL Lac object (PKS 0521-36) has sig- 
nificant dilution in both bands reducing the measured EWs. Its 
large non-thermal dilution is expected in terms of powerful jet 
emission from the BL Lac object. None of the inactive elliptical 
galaxies show any dilution in the A'-band indicating that the EWs 
of RGs and ellipticals are not significantly affected by non-stellar 
dilution. Another estimate of dilution can be obtained from com- 
parison of nuclear and ofF-nuclear continuum shapes, in which 
case redder spectrum is expected in nucleus. The comparison 
was done between off-nuclear and nuclear spectra of RGs, but 
there was no difference between the spectra. This indicates that 
no significant amount of dust is present in the nuclear region of 
RGs and, as a consequence, HII emission lines are not seen in 
the RG spectra. 

Table |4] gives the average values of ofF-nuclear and nuclear 
EWs of Sil, CO(1.62), Nal, Cal and CO(2.29) absorption Hues 
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for RGs and the corresponding values for inactive ellipticals. 
While the off-nuclear EWs of RGs are generally in agreement 
with the nuclear EWs of RGs and those of ellipticals, the off- 
nuclear EW(CO 2.29) of RGs is larger than that in the nuclear 
region of RGs and in elliptical galaxies. As was previously men- 
tioned, it is unlikely that this difference is caused by non-stellar 
dilution. The EW(CO 2.29) is sensitive to the luminosity class 
of stars being smaller in dwarfs than in giant stars. The differ- 
ence of EW(CO 2.29) between RGs and inactive ellipticals may 
indicate a different stellar populations. 

Figure |5] shows the EW(NaI+CaI) as a function of EW(CO 
2.29) for off-nuclear and nuclear RGs, and inactive elliptical 
galaxies together with stars of different spectral type. Caution 
should be exercised when comparing stellar and galactic data, as 
galaxies are composite stellar systems and their spectra show the 
integrated contributions from stars spanning a range of proper- 
ties. Still, given that the NIR spectral region is dominated by 
evolved (RGB) stars, this issue is not as critical as at visible 
wavelengths, where stars contribute over a much larger range 
of evoluti onary states (both RGB and main sequence stars e.g. 
iMarastonl ([2005)). All three luminosity classes of stars have 
a relatively tight linear correlation between EW(Nal-HCal) and 
EW(CO 2.29). Most of the off-nuclear RGs are in good agree- 
ment with the relation of giant (III) stars, while inactive ellipti- 
cals deviate more strongly from that relation. On the other hand, 
theplotof EW(CO 2.29) vs. EW(CO 1.62)/EW(Si 1.59) (Fig.© 
indicates that the stellar populations are mostly dominated by 
supergiants (I). However, it is possible that the CO(1.62) line is 
underestimated or the Sil line overestimated decreasing the line 
ratio of EW(SiI)/EW(CO 1.62) and thus, in both cases, shifting 
the real data points to the left in Fig.|6l It is reasonable that Sil 
line is slightly overestimated, because the absorption feature is 
located in the //-band which has numerous absorption lines and, 
in medium resolution spectra, the lines are blended. In conclu- 
sion, the majority of RGs are consistent with the relation of giant 
stars. Only a few RGs are readily in agreement with ellipticals, 
while the majority are displaced toward giants and supergiants 
indicating the difference in their stellar populations. 

Almost all RGs and inactive elhpticals have EW(CO 2.29)> 
15 A which is larger than the EW of K giants (from lOA to 
15A) or main sequence stars. In fact, the only stars with EW 
of CO(2.29) higher than 15A are M giants and K and M super- 
giants. The stellar population of these galaxies must therefore 
be dominated by very cool stars, probably M giants, to match 
the measured EWs. Note, however, that since M stars are much 
brighter in the /T-band than K stars, this does not imply that M 
giants are the only significant luminosity-weighted population. 

Star formation in the vicinity of nucleus can be activated by 
nuclear activity. Indeed, there is evidence for ongoing circum- 
nuclear sta r formation in the ce ntral region of the BL Lac PKS 
2005-489 dBressan et al.l l2006h . Different star formation rates 
within a galaxy might set up a radial stellar population gradient 
which can be seen as a difference in line strengths between nu- 
clear and off-nuclear spectra. However, non-stellar dilution must 
be taken into account that weakens the absorption lines of stellar 
population in the nuclear region. Only 3C 29 shows a systematic 
variation in the strength of the lines between off-nuclear and nu- 
clear spectra. The variation of the line strengths of 3C 29 can be 
due to a stellar population gradient but, however, deeper and bet- 
ter spatial resolution spectra are needed to confirm the possible 
gradient. 

EW(Cal) and EW(CO 2.29) of off-nuclear RGs and ellip- 
tical galaxies are plotted against EW(Nal) in Fig. I?] They are 



Table 4. Average values of EWs for off-nuclear and nuclear ra- 
dio galaxies and inactive elliptical galaxies compared to EWs of 
supergiants (SG), giants (G) and dwarfs (DW). 



Sil 



CO(6-3) Nal 



Cal 



CO(2-0) 



RG:off-nuc 3.7 ± 1.2 4.3: 

RG: nuc 2.6 ± 0.7 4.7 : 

Ellipticals: 2.9+ 1.4 4.9: 

Supergiants (I): 4.3±0.6 4.6h 

Giants (III): 3.8±04 5.1: 

Dwarfs (V): 3.2±0.8 1.8: 



tO.7 4.2 ±0.5 3.3: 

t 1.1 4.0 ±0.9 4.3: 

t 1.0 4.7 ±0.7 3.2: 
2.7 - 

t 1.9 2.8 ±0.5 2.9: 

tO.6 2.9 ± 1.9 2.6: 



0.3 20.3 ±44 
2.1 16.2 ±24 
1.0 16.6 ±2.3 
17.6±8.0 
0.8 164 ±4.7 
1.4 5.3 ±3.1 




PKS 0521-36 



0.2 
log EW(CO 1 .62)/EW(Si 1.59) 



Fig. 3. EW(CO 1.62) as a function of log[EW(CO 1.62)/(EW(Si 
1.59)] for off-nuclear spectra of radio galaxies (filled circles), 
nuclear spectra of radio galaxies (filled squares) and spectra of 
inactive elliptical galaxies (open triangles). Different luminosity 
types of stars are shown as symbols: supergiants I (open circles), 
giants III (asterisks ) from [R amirez e t alJ (II997) and dwarfs V 
(open squares) fromlAli et al. ( 1995). The lines enclose the area 
occupied by stars with no dilution and the arrow gives the direc- 
tion of the effects of dilution. 



compared with the relation of solar metallici ty cluster stars an d 
purely old population galaxies adapted from lSilva et al.l (l2008h . 
In the plot of Cal vs. Nal, majority of the ellipticals occupy the 
locus which is consistent with the relation of purely old galax- 
ies, whereas all RGs have larger Cal indices. Furthermore, in 
both EW(Cal) and EW(CO 2.29) plots, majority of the RGs and 
ellipticals also have larger Nal indices than those of solar metal- 
licit y cluster s tars, b eing consistent with the results obtained 
by .Silva et al.l (|2008|) for early-type galaxies of Fornax cluster 
Unlike found by Cesetti et al.l (l2008l) . our elliptical galaxy sam- 
ple do not show a significant correlation between EW(Cal) and 
EW(Nal). In both plots, RGs deviates from the ellipticals sug- 
gesting that elliptical galaxies with AGN are in a different stage 
in their evolution than inactive counterparts. 
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Fig. 5. The sum of the EWs of Nal and Cal as a function of 
EW(CO 2.29). The lines show the loci of different luminosity 
types of stars. The meaning of the symbols is the same as in 
Fig. [3] 



3.3. Age estimates for stellar population of radio galaxies 

It is well known that the majority of inactive elliptical galax- 
ies have not experienced recent significant star formation, but 
are dominated by passively evolved relatively old population 
of lat e-type giant RGB s t ars. However, many r ecent imaging 
(e.g., iGovoni et al.l l2000t iHy vonen et al.l l2007al) and spectro- 



scopic studies (e.g.. 'Raima nn et ani2005h have shown that, un- 
like inactive ellipticals, RGs have experienced a recent star for- 
mation episode and, thus, bluer host ga laxy colors are due to 
younger stellar population component. iRaimann et al.l (|2005|) 
obtained spectra for 24 RGs at optical wavelengths, and they 
found a systematic difference in stellar population between RGs 
and inactive ellipticals, in the sense that RGs have a larger con- 
tribution from an intermediate age population. 

Although stellar pop ulation synthesis models (e.g. , 
iFioc & Rocca-Volmerangg Il997t iBruzual & CharlotI l2003l) 
at NIR wavelengths are not as well developed as at visible 
wavelengths, they can be used to provide estimates for th e 
ages of the stellar population of galaxies dRiffel et al.l l2008l) . 
Theoretical spectral fitting is not currently possible at NIR, but 
synthesis models can be used to derive the time evolution of the 
EWs of diagnostic stellar lines parametrized by IMF and the 
metallicity of stellar population. 

To calculate the e volution of EW fo r the NIR absorption lines 
we used the model of Maraston (2005) with Salpeter IMF, solar 
metallicity and SSP model. That model was preferred over the 
continuous star formation model because the latter predicts Br- 
■y lines which are n ot seen in our RG spectra and, furtherm ore, 
in previous studies (lKennicuttlll998l:lBendo & Josephll2004 the 
observed spectra were found to be more consistent with the SSP 
model. Metallicity and age of stellar populations are degener- 
ate, in the sense that increasing metallicity generates older stel- 
lar populations. However, very low metallicities of RGs and el- 
lipticals are very unlikely because the observed line indices are 
significantly greater than those predicted by the low metallity 
models. 

The measured EWs of each absorption line of off-nuclear 
RGs and inactive elliptical galaxies were compared with corre- 
sponding EW vs. time theoretical curve computed by the synthe- 
sis model, allowing to estimate the ages of the stellar population. 
Note that before this comparison, the observed spectra (~30 A) 



8 



Hyvonen et al.: The stellar population of low redshift radio galaxies 



A 

A A 



A 

I I I I I 

3 4 

Nal 




pernova should be seen. The lack of these lines in the observed 
spectra implies either weak or nonexistent supernova activity 
and, thus, nonexistent young population. For these reasons, al- 
though both are formally consistent with the observations, the 
intermediate age population is preferred over the young popula- 
tion in Table |5] 

This result is in good agreement wit h the ages of n earby 
RGs derived from the optical spectra ( Raimann et al.ll2005 ). and 
suggests the widespread existence of a two-component stellar 
population containing both an old dominating component and 
an intermediate age component, in both RGs and inactive ellip- 
ticals. Note t hat a similar two- component model was recently 
suggested by iSilva et al.l (l2008l) for three early - type g alaxies of 
Fornax cluster. The sample of iRaimann et al.l (l2005b includes 
one common object to us, namely NGC 612. For that object, 
they estimated the stellar population to be ~ 2.4 Gyr of age that 
is consistent with our result (Table [S]). These ages are, on the 
other hand, only in reasonable agreement with those based on 
fitting single-age, single-metallicity stellar populatio ns to optical 
and NIR colours. For example, using this method, James et al] 
(2006) derived luminosity-weighted mean ages between 3 and 
14 Gyr, with a mean of ~8 Gyr, for field and Virgo cluster ellip- 
tical galaxies. However, we note that near-infrared is much more 
sensitive spectral region to trace evolved RGB stars than optical, 
in the sense that even a relatively small RGB fraction in a stel- 
lar population generates strong contribution to the strength of 
near-infrared absorption features. This causes a systematic ef- 
fect between the age estimates from optical and near-infrared 
methods. Although RGs and inactive ellipticals occupy different 
loci in the index vs. index plots of Fig. |2l reflecting differencies 
in their stellar content, the best-fit age estimates of both samples 
are in reasonable agreement. Most likely, this indicates different 
metallicities in RGs and inactive ellipticals. 

RGs appear to be slig htly older than late- type starburst galax- 
ies (~1 Gyr) observed bv lRiffel et al.1 (|2008|) in NIR. This might 
indicate that RGs have experienced nuclear activity and en- 
hanced star formation earlier than starburst galaxies, and in those 
RGs star formation have finished possibly due to AGN feedback 
phenomenon. 



Fig. 7. EWs of Cal and CO(2.29) against EW(Nal). Symbols are 
the same as in Fig. [3] The dotted and long-dashed lines represent 
the relations of the solar m etallicity cluster star fits and purely 
old population (ISilva et al.ll2008) . respectively. 



were smoothed to the same resolution of the MarastonOS model 
(100 A). Theoretical curves were calculated using exactly the 
same parameters as th ose in a r ecent NIR spectroscopic study 
of starburst galaxies bv lRifi'el eT al. (2008), and we refer to their 
Fig. 4 for the behavior of the absorption lines as a function of 
age. Fitting was performed for /T-band Cal feature and CO(2.29) 
bandhead, and the results of the fitting are given in Table|5] Note 
that the theoretical EW values are not unique with evolutionary 
time, but the same EW can have two to four diff'erent time values 
between very young (~10 Myr) and intermediate age (~2 Gyr) 
stellar population fits. EWs of Cal and CO(2.29) absorption lines 
of RGs are formally consistent with the fits of a young (~400 
Myr) and intermediate age (~2 Gyr) stellar population. However, 
the young population is unlikely because massive stars passing 
the red giant phase will explode as supernovae, and emission 
lines of [FeII](1.644yum) and H2l-0S(l)(2.212;um) from the su- 



4. Conclusions 

We have observed the first NIR HK-hand spectra of a sizable 
sample of low redshift RGs to characterize the EWs of their di- 
agnostic NIR absorption line features and, in addition, based on 
the evolution models of the diagnostic features trace the inter- 
mediate age stellar population component. RG sample was com- 
pared with nearly similar size inactive elliptical galaxy sample 
obtained with the same telescope, instrument and wavelength 
range. For both samples spectral indices of several absorption 
features were measured in both bands, namely Sil (1.589//m), 
CO(6-3) (1.619//m) in the //-band, and Nal (2.207//m), Cal 
(2.263yum), CO(2-0) (>2.29yum) in the /T-band. These fines are 
sensitive to the stellar temperature and/or luminosity so they can 
be used as a tool for identifying different spectral and lumi- 
nosity classes of stellar population. NIR spectroscopy is pow- 
erful to study intermediate age stellar populations having advan- 
tages over optical spectroscopy to trace TP-AGB stars. Although 
currently stellar population synthesis models are not fully self- 
consistent for the interpretation of NIR spectra of galaxies, they 
still can be used to estimate the age of stellar population by com- 
paring the measured EWs of diagnostic stellar lines to the EWs 
calculated by theoretical evolution curve. 
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Table 5. Age estimates in Myr for the stellar population com- 
ponent of each off-nuclear radio galaxy and inactive elliptical 
galaxy together with the stellar population age of the combined 
RG and inactive elliptical samples." 



Target 


Cal 




CO 


(2.29) 


Young 


Interm. 


Radio Galaxies: 














PKS 0023-33 


350 


1580 


280 


2000 


315 


1800 


3C 29 


340 


1800 


680 


1260 


510 


1500 


NGC 612 


320 


2000 


630 


1580 


480 


1800 


ESQ 552-G14 


320 


2000 


320 


2630 


320 


2300 


PKS 0521-36 


280 


2240 






280 


2200 


ESQ 528-G36 


310 


2000 


700 


1300 


510 


1650 


NGC 6998 


350 


1600 






350 


1600 


PKS 2158-380 


330 


1800 


700 


1300 


500 


1500 


Combined RG 


350 


2000 


660 


1400 


500 


1700 


Ellipticals: 














NGC 2640 


320 


2000 


700 


1300 


510 


1650 


NGC 4478 


360 


1600 


270 




320 


1600 


NGC 4472 


280 


2140 


400 


2000 


340 


2100 


NGC 4546 


260 


2240 


380 


2040 


320 


2100 


NGC 4564 


280 


2140 


400 


2000 


340 


2100 


NGC 4621 


630 


1100 


370 


2400 


500 


1750 


NGC 4649 


280 


2100 


370 


2100 


330 


2100 


NGC 4697 


340 


1800 


360 


2500 


350 


2150 


NGC 5576 


260 


2300 


320 


2800 


290 


2550 


Combined E 


340 


1900 


350 


2200 


350 


2050 



" Column (1) gives the name of target; (2) and (3) best-fit young and 
intermediate ages based on the EWs of Cal; (4) and (5) ages based 
on the EWs of CO(2.29), and (6) and (7) the average age of the 
young and intermediate age stellar population component. 



Plotting EW(Nal-i-Cal) vs.EW(CO 2.29) RGs are consistent 
with the relation of giants stars, whereas inactive ellipticals de- 
viates more strongly from the relation toward main sequence 
dwarf stars, and on average, EW(C02.29) of RGs is somewhat 
greater than that of inactive ellipticals. Most likely, EWs are 
not significantly affected by non-stellar dilution, and thus it in- 
dicates that elliptical galaxies containing AGN are in a differ- 
ent stage in their evolution than their inactive elliptical counter- 
parts. Our result supports previous optical spectroscopic studies 
of RGs which have shown evidence from intermediate age stel- 
lar population component, indicating a link between a nuclear 
activity and star formation in the host galaxy. Furthermore, RGs 
appear to contain a slightly older stellar population component 
than those of late-type starburst galaxies suggesting that RGs 
have experienced enhanced star formation earlier than starburst 
galaxies. There are no H2l-0S(1) and [Fell] emission lines in 
the RG spectra indicating either only very small or non-existing 
contribution from very young stellar population. 

To have a better understanding of the stellar content of RGs, 
larger NIR spectroscopic sample is needed with higher spectral 
and spatial resolution. This gives the possibility to obtain veloc- 
ity dispersion of the galaxies which, together with luminosity, 
can be used to calculate the M/L ratio of galaxies that is essen- 
tial information for tracing the young population. Furthermore, 
the correlation between velocity dispersion and K-band Nal in- 
dex can be used as a young population i ndicat or for early-type 
galaxies as was suggested bv lSilva et al.l (l2008l) . To characterize 
stellar population gradients from the central region, where AGN 
and star formation may co-exist, to the outer regions of galaxies 
higher spectral and spatial resolution are needed. Furthermore, 



to better interpret the composite NIR spectra of galaxies, more 
extended stellar population synthesis models should be devel- 
oped. To assert the unified model of radio-loud AGNs, where 
RGs are the parent population of BL Lac objects, larger samples 
of NIR spectra of BL Lac objects are required together with high 
resolution optical spectroscopy. 
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